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Ch Abstract 

. rH The CMS silicon tracker, consisting of 1440 silicon pixel and 15 148 silicon strip de- 

tector modules, has been aligned using more than three million cosmic ray charged 
particles, with additional information from optical surveys. The positions of the mod- 
ules were determined with respect to cosmic ray trajectories to an average precision of 
3-4 microns RMS in the barrel and 3-14 microns RMS in the endcap in the most sen- 
sitive coordinate. The results have been validated by several studies, including laser 
beam cross-checks, track fit self-consistency, track residuals in overlapping module 
regions, and track parameter resolution, and are compared with predictions obtained 
from simulation. Correlated systematic effects have been investigated. The track pa- 
rameter resolutions obtained with this alignment are close to the design performance. 
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1 Introduction 

The main goal of the CMS experiment [1] is to explore physics at the TeV energy scale exploit- 
ing the proton-proton collisions delivered by the Large Hadron Collider (LHC) [2 J. The CMS 
silicon tracking detector (tracker) ||3j [4J consists of 1440 silicon pixel and 15 148 silicon strip de- 
tector modules and is located, together with the electromagnetic and hadronic calorimeters, in- 
side a superconducting solenoidal magnet operating at a field of 3.8 T. Outside of the solenoid, 
the muon system is used both for triggering on muons and reconstructing their trajectories. 

The aim of the tracker is to measure the trajectories of charged particles (tracks) with excellent 
momentum, angle, and position resolution, and with high pattern recognition efficiency 0. 
Precise determination of the position of all silicon modules is a challenging task and one of the 
critical aspects for achieving the design track parameter resolutions. In the context of align- 
ment, position is used throughout this paper to refer both to the location of the center point of 
active areas on modules and to the orientation of active areas. The hit position resolution is 
in the range 10 to 30 }im and therefore the alignment precision should be better than 10 }im 
to achieve optimal track parameter resolution. Simulation studies based on the design (ideal) 
tracker geometry imply an impact parameter resolution of about 15 and a transverse mo- 
mentum (pr) resolution of about 1.5% for 100 GeV/ c muons. 

The CMS collaboration conducted a month-long data-taking exercise known as the Cosmic Run 
At Four Tesla (CRAFT) during October-November 2008 with the goal of commissioning the ex- 
periment for extended data-taking (6j. With all installed detector systems participating, CMS 
recorded 270 million cosmic-ray-triggered events with the solenoid at its nominal axial field 
strength of 3.8 T. However, only a few percent of those events had cosmic rays traversing the 
tracker volume. Prior to CRAFT and during the final installation phase of the experiment, a se- 
ries of commissioning exercises to record cosmic ray events took place with the solenoid turned 
off from May to September 2008. The tracker was included only at the later stages of those runs 
and a smaller number of cosmic ray events with the tracker participation were recorded com- 
pared with CRAFT. The operating temperature of the tracker during the CRAFT data- taking 
period was stable at around room temperature. Excellent performance of the tracking system 
has been achieved with both the silicon strip [7j and silicon pixel [8] components. 

In this paper, the measurement of detector module positions for the full tracker from alignment 
analysis with cosmic ray particles, optical survey information, and the Laser Alignment Sys- 
tem (LAS) is presented. This work builds on the experience with the alignment of about 15% 
of the CMS silicon strip tracker during stand-alone commissioning [9j. Significant improve- 
ments have been introduced. Inclusion of the pixel detector and of the full strip tracker made 
the first alignment of the whole CMS tracker possible. The magnetic field allowed estimates 
of the multiple scattering effects on a track-by-track basis. The data taken with the solenoid 
turned off were used in this analysis to cross-check the field-on data in a small number of 
cases. Significantly more tracks in a single stable configuration improved statistical precision. 
Two complementary statistical alignment algorithms, a global and a local iterative method, 
were applied to the data. Although both algorithms had been used previously during stand- 
alone commissioning, further improvements have been introduced, such as inclusion of survey 
data together with tracks. Results from the alignment studies have been compared extensively 
with detailed Monte Carlo (MC) simulations of the detector. 

This paper is structured as follows: in Section |2j the tracker layout and coordinate systems are 
introduced. Section [3] summarizes the optical survey measurements. In Section]!} the track re- 
construction and alignment algorithms are described and their application to data is explained. 
Results are presented in Section |5j 
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2 Tracker Layout and the Coordinate Systems 
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Figure 1: A quarter of the CMS silicon tracker in an rz view. Single-sided silicon strip module 
positions are indicated as solid light (purple) lines, double-sided strip modules as open (blue) 
lines, and pixel modules as solid dark (blue) lines. Also shown are the paths of the laser rays 
(R), the beam splitters (B), and the alignment tubes (A) of the Laser Alignment System. 



2 Tracker Layout and the Coordinate Systems 

CMS uses a right-handed coordinate system, with the origin at the nominal collision point, 
the x-axis pointing to the center of the LHC, the y-axis pointing up (perpendicular to the LHC 
plane), and the z-axis along the anticlockwise-beam direction. The polar angle (6) is measured 
from the positive z-axis and the azimuthal angle ((p) is measured from the positive x-axis in the 
x-y plane, whereas the radius (r) denotes the distance from the z-axis. 

The CMS tracker consists of a silicon pixel detector and a silicon strip detector (Fig. [IJ. The 
silicon pixel detector is composed of two sub-detectors, the barrel (BPIX) and the two endcaps 
in the forward regions (FPIX). The pixel modules provide two-dimensional measurements of 
the hit position in the module planes, which effectively translate into three-dimensional mea- 
surements in space. The silicon strip detector is composed of four sub-detectors: the Tracker 
Inner and Outer Barrels (TIB and TOB), the Tracker Inner Disks (TID), and the Tracker End- 
caps (TEC). All sub-detectors are concentrically arranged around the nominal beam axis. The 
two inner layers of both the TIB and TOB, the two inner rings of the TID, and the first, sec- 
ond, and fifth rings of the TEC are equipped with double-sided modules; all other positions 
have single-sided modules. Single-sided modules provide rty measurements in the barrel and 
(p measurements in the endcaps. Double-sided modules are made of a pair of single-sided strip 
modules, one rep and one stereo module in the barrel, and one (p and one stereo module in the 
endcaps, precisely mounted back-to-back, with the stereo module sensors tilted by 100 mrad. 
For simplicity, we refer to both rep and (p modules as rep in the rest of the paper. 

In the barrel region, modules are arranged in linear structures parallel to the z-axis, such as 
ladders in the BPIX, strings in the TIB, and rods in the TOB. The endcaps are composed of 
disks, which in turn contain wedge-shaped structures covering a narrow cp region, such as 
blades in the FPIX and petals in the TEC. The BPIX and the TIB are composed of two half- 
barrel structures, separated along the x = plane for the BPIX and the z = plane for the 
TIB. 

A local right-handed coordinate system is defined for each module with the origin at the ge- 
ometric center of the active area of the module. As illustrated in Fig. |2j the w-axis is defined 
along the more precisely measured coordinate of the module (typically along the azimuthal 
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direction in the global system), the c-axis orthogonal to the w-axis and in the module plane, 
pointing away from the readout electronics, and the w-axis normal to the module plane. When 
double-sided modules are considered as a single entity, the coordinate system is referenced to 
the rep module. For the pixel system, u is always orthogonal to the magnetic field, that is in 
global rep direction in the BPIX and in the radial direction in the FPIX. The v coordinate is per- 
pendicular to u in the sensor plane, that is along global z in the BPIX and at a small angle to the 
global r(p direction in the FPIX. The angles a, /S, and 7 indicate rotations about the u, v, and w 
axes, respectively. 

In addition, local u' and v' coordinates are defined such that they are parallel to u and v, but the 
direction is always chosen to be in positive (p, z, or r directions, irrespective of the orientation of 
the local coordinate system. For the TID and TEC wedge-shaped sensors, where the topology 
of the strips is radial, the u'- and o'-axes change direction across the sensor such that v' is always 
directed along the strips and therefore u' corresponds to the global r^-coordinate. 




Figure 2: Illustration of the module local coordinates u, v, w and the corresponding rotations 
a, f>, 7 for a single-sided strip module. 



3 Optical Survey Measurements 

Optical surveys taken during module construction and integration provide initial alignment 
parameters for many of the modules. Additionally, the survey information was used as a con- 



straint in the alignment procedure, as will be described in Section 4.1 



Coordinate Measuring Machine (CMM) data and photogrammetry have been used for the op- 
tical survey of the tracker components. While the former were used for measurements of the 
active elements, the latter were used for the alignment of larger structures. For the inner strip 
detectors (TIB and TID), both module-level and high-level structure information were used in 
the analysis. For the outer strip detectors (TOB and TEC), module-level survey data were used 
only as a conformance criterion during construction, while survey information from high-level 
structures was used in the analysis. More information on the strip detector surveys can be 
foundinRef.il. 

Detailed optical surveys of the pixel endcap detectors were performed as part of the construc- 
tion process. First, module positions were measured within a panel, which contains three or 
four modules. Then the positions of modules were also measured on a half-disk, where 12 
panels are placed on each side. Eight fiducial points were measured on the active area of each 
module, allowing a position measurement precision of about one micron. Some of the module 
fiducial points were partially obscured by other modules when mounted on a half-disk and 
thus could not be measured. Good reproducibility of module measurements within a panel 
was found in the two configurations, allowing an estimate of survey precision. Measurements 
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4 Track-based Alignment 



of the two sides of the half-disk were connected through a touch-probe, which related posi- 
tions of the touch-probe targets previously measured with CMMs. This also related measure- 
ments from the front and the back panels of blades. Finally, half-disks were placed in the half- 
cylinders, which were then inserted in the pixel detector volume. Half-disk and half-cylinder 
positions were measured with photogrammetry and were related to the positions of the active 
elements through photo targets, which had been previously measured with CMMs. 

In the barrel pixel detector, only two-dimensional measurements of the module positions within 
a ladder were performed. High-resolution digital images of four nearby fiducial points in 
neighboring modules were taken. A precision of 2 /mi was achieved for longitudinal posi- 
tions of the modules within a ladder, allowing, in particular, a tight constraint on the z-scale. 
The relative position of nearby ladders in the half-shells were also surveyed with similar two- 
dimensional measurements, but were not used in this analysis. 



4 Track-based Alignment 

The goal of the track-based alignment procedures is to determine the module positions from a 
large sample of reconstructed charged particle trajectories. Each trajectory is built from charge 
depositions ("hits") on individual detectors, assuming a piece-wise helical track model, incor- 
porating effects from multiple scattering and energy loss. The "Combinatorial Track Finder " 
(CTF) track algorithm [7J was used to reconstruct the cosmic muon trajectories. The result is 
a set of track parameters. Five track parameters describe the helical trajectory of a track at the 
point closest to the nominal interaction point: distance of closest approach in the transverse 
d xy and longitudinal d z directions, track azimuthal angle cp, track polar angle 9, and transverse 
momentum pj. 
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Figure 3: Momentum (left), azimuthal angle (middle), and polar angle (right) spectra of cosmic 
muons reconstructed in the CMS tracker volume based on the selection criteria described in the 
text. The solid (red) circles represent the cosmic ray data whereas the open (blue) circles come 
from a MC simulation. The vertical cosmic ray tracks correspond approximately to cp = 6 = 
90°. 



The momentum and angular spectra of cosmic muons reconstructed in the CMS tracker volume 
are shown in Fig. |3j where good agreement with MC simulation [10J is observed. These track 
parameters are defined at the point of closest approach to the CMS nominal interaction point 
and the absolute value of the (p angle is shown, considering all particles traversing CMS from 
top to bottom. The aligned tracker geometry was used in track reconstruction. The asymmetry 
about cp = 90° is due to the excess of positive over negative cosmic ray particles entering CMS 
and the bending of charged tracks in the magnetic field of the solenoid. Positively charged 
cosmic rays tend to peak just above 90° whereas negatively charged cosmic rays tend to peak 
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just below 90° . The asymmetry about 6 = 90° is mainly due to the asymmetric location of the 
cavern shaft. Tracks were selected according to the following criteria. Each track was required 
to consist of at least eight hits with a signal-to-noise ratio larger than 12 in the silicon strip 
modules or a probability of the pixel hit matching the template shape [11] of at least 0.001 (0.01) 
in the local u (v) direction. Hits were also rejected if the track angle relative to the local uv plane 
was less than 20°. In addition, two hits were required to be on either pixel or double-sided strip 
modules, allowing a precise measurement of the polar angle 8. The particle momentum had to 
be greater than 4GeV/c, and the x 2 P er degree of freedom of the track fit, x 2 /ndf, had to be 
less than 6.0. Outlying hits with large normalized residuals were removed. In total, about 3.2 
million tracks were selected for alignment, out of which about 110000 had at least one pixel 
hit. 

The track reconstruction used the Alignment Position Errors (APE), the estimated uncertainty 
on the module position in the three global coordinates, which were added in quadrature to 
the hit errors during the pattern recognition and track fitting procedure. This allowed efficient 
track reconstruction in the presence of misalignment and a correct pull distribution of track 



parameters. The determination of APEs is described in Section 5.1 The APE values used for 
the initial track reconstruction before alignment analysis were set to large values as they had to 
account for possible large displacements of the entire sub-detectors while still guaranteeing an 
efficient track-hit association. 



4.1 Alignment algorithms 

The alignment is an optimization problem that can be formulated in the context of linear least 
squares. Module position corrections ("alignment parameters") p are determined by minimiz- 
ing an objective function 

tracks hits 

which can be expressed as the sum over all hits i on all tracks j and track parameters qy, as- 
suming negligible correlations between hits. Track residuals tu = m» — tu (p, qy) are defined as 
the difference between the measured hit position m« and the trajectory impact point f„-, both 
of which are given in the module local coordinate system. The residuals are either one- or two- 
dimensional vectors, depending on the type of module, and V„- is either the squared error or 
the covariance matrix in case of one- or two-dimensional residuals, respectively. 

Two statistical methods were employed to solve the alignment problem. Both of them were 
previously applied to the CMS silicon strip tracker alignment during stand-alone commission- 
ing 0. The global alignment algorithm ("Millepede II") lTl2l| minimizes the % 2 function in 
Eq. ^ by taking into account track and alignment parameters simultaneously. This algorithm 
has also been previously used in simulation studies of the CMS tracker |T3| . Since only the 
n alignment parameters p are of interest, the problem is reduced to the solution of a matrix 
equation of size n. The value of n is of the order of 10 5 for six degrees of freedom of 16588 
modules. The covariance matrix V,-; is approximated by a diagonal matrix with uncertainties 
Cij for uncorrected hit measurements m„- of the track Given reasonable starting values, po 
and qyo, the track model prediction fy(p, qy) can be linearized. Since angular corrections are 
small, the linearized problem is a good approximation for alignment. 

The local iterative algorithm ("Hits and Impact Points") |fl4l approximates Eq. jl]) by assuming 
no track parameter q dependence and therefore ignores correlations between alignment pa- 
rameters for different modules in one iteration. The trajectory impact point f a y is recalculated 
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4 Track-based Alignment 



for each hit, removing the hit under consideration from the track fit. The track parameters and 
correlations between different modules are taken into account through iterations of the mini- 
mization procedure and refitting the tracks with new alignment constants after each iteration. 
The number of iterations is determined empirically. This approach allows a simplification of 
the formalism and leads to a solution of the six-dimensional matrix equation for the six align- 
ment parameters of each module in each iteration. 

The local iterative algorithm permits the inclusion of survey measurements in the formalism of 
Eq. jlj, as described in Ref. |Tl5l . This leads to an additional term in the objective function to be 
minimized independently for each module in a given iteration: 

hits survey 

xKp) = E r f(p) v rVp) + E r ',(p) v ;/ r *,(p)' ( 2 ) 

where the original formulation has been modified to include survey information from the k hi- 
erarchial levels for each sub-detector. The track residuals f;(p) do not have explicit dependence 
on track parameters and enter the sum over hits in a given module m. The six-dimensional sur- 
vey residuals r # i are defined as the difference between the reference and the current sensor 
position. The survey measurement covariance matrix V*y reflects both the survey precision 
and additional uncertainties due to changes in the detector. These errors can be configured 
differently for different hierarchy levels and for the degrees of freedom that should be stable, 
such as the longitudinal direction in a barrel ladder, and those which may change more with 
time. 

The local iterative method uses the full implementation of the Kalman filter track reconstruc- 
tion algorithm adopted in CMS 0- It requires a large number of iterations and large computing 
resources to refit the tracks in each iteration. The global method, instead, allows the determi- 
nation of alignment parameters, properly accounting for the correlations among them, in a 
single step. However, the global method does not take into account the effects of material in 
the tracker and assumes a simple helical trajectory for charged particles. The method also re- 
quires a large amount of computer memory and the application of advanced techniques |12| 
for solving Eq. qTp . Each of the two alignment algorithms was used to obtain module positions 
independently and a comparison of results between the two complementary approaches was 
part of the validation of the tracker alignment procedure. After verifying that the two methods 
yielded consistent results, the final results were obtained by applying the two algorithms in 
sequence in order to take advantage of their complementary strengths. 

4.2 Alignment strategy 

Single-sided silicon strip modules provide only a one-dimensional measurement in the module 
plane, along the local u-coordinate, which effectively translates into a two-dimensional mea- 
surement in space. The ^-coordinate is only known to be within the module boundaries, with 
precision not sufficient for track reconstruction requirements. However, the information from 
the rep and stereo modules in a double-sided module is combined into a two-dimensional mea- 
surement in the combined module plane in both u and v for the pattern recognition phase. 

Due to the 100 mrad stereo angle between the rep and stereo modules, though, a small displace- 
ment in u is equivalent to a ten times larger displacement in v. Given comparable mounting 
precision of modules in u and v, a smaller u displacement is much more likely than a larger v 
displacement. In fact, it was found from the data after several attempts at alignment of double- 
sided modules in v, that parameters much larger than the known assembly accuracy were ob- 
tained. Whenever possible, the two single-sided components of a double-sided module were 



4.3 Alignment with the global method 



7 



aligned separately, but only in the most precise coordinates. The relative alignment of the two 
sides was found to be consistent with assembly tolerance and at the same time improved track 
residuals significantly. This fact and the hierarchical structure of the tracker led to multi-step 
strategies for both alignment methods. 

4.3 Alignment with the global method 

The first step of the global method aligns the highest-level structures (half-barrels, endcaps) 
with all six degrees of freedom together with all module units, including rep and stereo strip 
modules in a double-sided module, with the most sensitive degrees of freedom each (u, w, 7, 
and for pixel modules also v). The detector design geometry was chosen as a starting point. A 
limitation of 46 340 parameters in the program used led to selection criteria for modules to be 
aligned of more than 25 hits per pixel module and more than 425 hits per strip module. This 
resulted in the alignment of about 90% of all modules. 

In the second step all modules (all double-sided or single-sided strip modules with more than 
150 hits and all pixel modules with more than 25 hits) were aligned in the TIB in u, w, a, f>, 7; 
in the pixel system in u, v, w, 7; and in u, w, 7 elsewhere. Compared to the previous step, this 
recovered the remaining 10% of the modules and allowed more degrees of freedom for the TIB, 
which had larger assembly tolerance, but did not allow independent alignment of the rep and 
stereo modules within double-sided combinations. No alignment of a and /3 was performed 
in the TOB due to its higher mounting accuracy, and in the TID and TEC, since they were less 
often traversed by cosmic ray tracks. 

The third step was designed to recover lost correlations between the first two steps and had the 
same configuration as the first step, but the minimum number of hits in the strip modules was 
increased to 450 with respect to 425 used in the first step. 

4.4 Alignment with the local method 

Module positions as determined from optical surveys were used as the starting point in the lo- 
cal method. The difference of the starting geometry to that of the global method was motivated 
by the use of survey information in the minimization procedure, as described in Section |4.1| 
However, for both methods the initial geometry did not affect the results significantly. The 
alignment with the local algorithm started with the large structures and then proceeded in or- 
der of increasing granularity down to the module units in order to speed up convergence. The 
small fraction (4%) of tracks passing through the pixel detector suggested a modified approach 
in which the silicon strip modules were first aligned without pixel information. The hits from 
the structures not yet aligned were excluded from the track fit. The APE values were set high 
at the beginning of the alignment process (several hundreds of microns depending on the sub- 
detector), and were progressively reduced to a few tens of microns in subsequent alignment 
steps. The values of the APEs were kept fixed within each alignment step. 

The alignment with the local method was carried out in six steps. In the first step (15 itera- 
tions), the half-barrels of BPIX, TIB, and TOB were aligned as rigid bodies with four degrees 
of freedom (u,p,w,7). In order to have a better sensitivity along the global z-direction, only 
tracks passing through the pixel barrel system were used. In the second step (15 iterations), 
the endcaps of the FPIX, TID, and TEC were aligned in the same four degrees of freedom. The 
third alignment step (30 iterations) moved the strip modules (treating double-sided modules 
as rigid bodies) in all six degrees of freedom. Data from the optical survey measurements (or 
design geometry if no survey information was available) were used according to Eq. The 
information on the module positions within a higher-hierarchy structure, for example within a 
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string of TIB, was used in the x 2 minimization. This proved to be useful for both limiting the 
movement of modules in poorly constrained degrees of freedom and aligning modules that 
had fewer than the required number of hits. In the fourth step (30 iterations), the strip module 
units (treating r(p and stereo strip modules in a double-sided module independently) with at 
least 50 hits were aligned in three degrees of freedom (u, w, 7). In step five (15 iterations), the 
ladders of the pixel detector were aligned in six degrees of freedom. Finally in the sixth step (15 
iterations), the pixel modules with at least 30 hits were aligned in all six degrees of freedom. 

After every step the alignment parameters were checked for convergence. Fewer than 20 mod- 
ules failed to converge since they were moved to implausible positions by the algorithm. The 
alignment parameters for these modules were left at the survey values. It was found that in 
some cases the basic procedure needed to be applied a second time in order to improve the 
convergence for the alignment of the large-scale structures. However, as will be discussed in 
the next section, the final alignment results were derived from a sequential approach involving 
both the global and the local method, and for this reason the local method was limited to a 
single set of steps. 

4.5 Alignment with the combined method 

The strength of the global method is solving effectively the global correlations, while for the 
local method it is that the same track fit is used as in the standard CMS reconstruction, and that 
survey information can easily be incorporated thus allowing for alignment with more degrees 
of freedom. In order to take advantage of both methods, the final alignment parameters were 



produced starting from the output of the global method analysis described in Section 4.3 then 



further aligning the tracker with the local method presented in Section 4.4 



The alignment strategy adopted by the local algorithm was modified to exploit the already 
good starting position of the modules provided by the global method. The APEs were kept 
fixed at 10 f/m in the barrels and 30 /mi in the endcaps, similar to those described in Section [BTT] 
In the first step of the local method (30 iterations), all strip modules (treating double-sided 
modules as rigid bodies) were aligned in six degrees of freedom using track and survey infor- 
mation. In the second step (20 iterations), the strip module units (treating r<p and stereo strip 
modules in a double-sided module independently) were aligned in three degrees of freedom 
(u,w,j). Pixel modules were not aligned in the first two steps, although pixel hits were included 
in the track fit. Finally, in the last step (20 iterations), the pixel modules were aligned in six de- 
grees of freedom. Modules for which the fit did not converge were left at the position obtained 
by the global algorithm. 



5 Results and Validation 

The alignment analysis was performed after the proper calibration of both the strip and pixel 
detectors 0IID- In particular, an appropriate Lorentz angle calibration was applied. A Lorentz 
angle miscalibration would result in a systematic shift of the effective charge collection posi- 
tion, in the same u-direction for all the modules with the same relative orientation between 
the solenoid and the drift field. Since this would appear as an effective shift of the module, it 
would not be evident if corrected for by the alignment. Therefore, for validation purpose only, 
alignment results were checked with the data collected both with T and with 3.8 T solenoid 
fields. The precision achieved in the alignment analysis is such that there should be clear evi- 
dence for a Lorentz angle miscalibration with an effective shift of only a few microns. In fact, 
good agreement with the Lorentz angle calibration [7] was found. 



5.1 Calibration of alignment position errors 
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Several approaches were employed to validate the alignment results. The low-level quanti- 
ties that were used in the x 2 minimization, such as residuals and the ^; 2 /ndf of the tracks, were 
monitored. Given the limited size of the alignment track sample, it was used both for alignment 
and validation. In addition, a high-level validation analysis was performed to monitor the track 
parameter resolution. Moreover, techniques for comparing positions of modules within differ- 



ently derived geometries, described in Sections 5.5 and 5.6 allowed a better understanding of 
the alignment performance. 

5.1 Calibration of alignment position errors 

In the track refit the APEs for each hit were assumed to be the same in the three spatial di- 
rections. For an alignment performed using cosmic ray tracks, the precision of the alignment 
varies within a given sub-detector because of the different illumination of modules due to their 
orientation relative to cosmic rays. The radius of the sphere representing the APE for each 
module, r sp here/ was therefore taken to be r sp h er e = r o/ v^entries/ where N en tries represents the 
number of hits in a module. The value of ro was chosen in order to have the Gaussian standard 
deviation of the distribution of the residuals normalized to their error approximately equal to 
one in the symmetric interval covering 95% of the distribution. In the TIB and TOB, a common 
yq value was defined for each layer and separately for the rep and stereo components of the 
double-sided layers. Elsewhere the value of ro was defined at the sub-detector level. The val- 
ues of r sp here were restricted to reasonable values, especially in the case of small or zero N en t r i es , 
where a precision compatible with survey and assembly data was used. APEs were calibrated 
using the set of alignment parameters obtained from the combined method. Examples of the 
distribution of the normalized residuals after the calibration of the APEs are shown in Fig. |4j 
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Figure 4: Distributions of normalized track residuals after the APE calibration procedure: BPIX 
(left), TIB (center), and TOB (right). Solid lines represent the results of Gaussian fits and the fit 
mean and sigma values are given within the plots. The alignment obtained from the combined 
method was used. 



5.2 Track fit quality and hit residuals 

All tracks used in the validation procedure were refitted with APEs corresponding to the align- 
ment parameters obtained with the combined method. The loose track selection described in 
Section [4] was applied. To avoid a bias in the determination of the hit residuals, the track pre- 
diction was calculated without information from the hit under consideration. The track ^ 2 /ndf 
distribution is shown in Fig. |5j The hit residuals in the u' and v' directions are shown in Fig. [6] 
The hit residual width is dominated by two effects other than alignment: track extrapolation 
uncertainties due to multiple scattering and hit position reconstruction uncertainties. Both of 
these effects are random, while misalignment leads to systematic shifts of the residuals. For this 
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Figure 5: Distributions of the ;f 2 /ndf of the tracks before alignment (dotted line) and after 
alignment with the local (dashed-dotted line), global (dashed line), and combined (solid line) 
methods. 



reason, the distribution of the median of the residuals (DMR) is taken as the most appropriate 
measure of alignment. Median distributions are shown in Fig. [7] and the corresponding RMS 
values of these distributions are given in Table [I] 

To check the statistical precision of track-based alignment a MC simulation was performed 
in which module positions from the combined method obtained with data were used as the 
starting geometry in the MC alignment procedure. This approach in MC effectively models 
the situation in data prior to and during the alignment. The resulting DMRs are also shown in 
Fig. [7] and the RMS values listed in Table [T] For comparison, the distributions obtained from 
the ideal MC simulation are presented in Fig.|7]as well. 

The middle columns of Table [T] reflect different stages of alignment: the initial situation with no 
alignment, alignment with tracks using the global and local methods, and finally with the com- 
bined method. Table[T]also includes results from MC simulations based on combined and ideal 
geometries. For the 2D pixel detectors both u' and v' coordinate measurements are quoted. 
Overall, there is significant improvement in the track reconstruction going from the geome- 
try without any alignment, to the alignment using tracks with the local and the global method, 
and finally to the combined result. With respect to cosmic ray trajectories, the module positions 
were determined to an average precision of 3-4 }im RMS in the barrel and 3-14 ^m RMS in the 
endcaps in the most sensitive coordinate. These values are in agreement with the expected 
statistical precision as determined using simulated events. They are also comparable to values 
obtained from a MC simulation based on the ideal detector geometry. 



5.2 Track fit quality and hit residuals 
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Figure 6: Track residuals, shown for BPIX (top left u' , top right v'), TIB (bottom left), and 
TOB (bottom right). The four lines correspond to positions before alignment (dotted lines) and 
after alignment with the global (dashed lines), local (dot-dashed lines), and combined methods 
(solid lines). 



Table 1: RMS of the distribution of the median of the residuals (DMR) in the u' and v' local 
coordinates for modules with more than 30 hits. The number of these modules compared to the 
total number of modules is stated in the last column. Four geometries are considered: those 
obtained with the three methods discussed in the text and the geometry before alignment. 
Results from simulations based on the combined alignment and ideal geometries are shown 
for comparison. 
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Figure 7: Distribution of the median of the residuals, for modules with more than 30 hits, 
shown for BPIX (top left u', top right v'), FPIX (second row left u' , second row right v 1 ), TIB 
(third row left), TOB (third row right), TID (bottom left), and TEC (bottom right). Shown are 
distributions before alignment (black dotted), after alignment with the combined method (red 
solid), combined method MC (green dashed), and ideal MC (blue dash-dotted). 



5.3 Residuals in overlapping module regions 
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5.3 Residuals in overlapping module regions 

A further method to monitor and validate the results of the alignment is to use the hits from 
tracks passing through regions where modules overlap within a layer of the tracker. This 
method, described in detail in Ref. |[16|, is also used to measure the hit resolution of the sen- 
sors [7J. In this method, the difference in residual values for the two measurements in the 
overlapping modules is compared, once the hits in the layer under consideration are removed 
from the track fit. The proximity of the hits reduces the amount of material between the two hits 
and minimizes the uncertainties on the predicted positions associated with track propagation. 
Deviations between the reconstructed hits and the predicted positions allow an assessment of 
the relative alignment between two adjacent modules. A non-zero mean of a Gaussian fit to the 
distribution of residual differences indicates a relative shift between the modules. A relative 
rotation between the modules can be seen in the dependence of the differences between the 
residuals and the track position predictions on the modules under study, where the slope indi- 
cates the magnitude of the residual rotation. This would also lead to an increase in the width of 
the Gaussian fitted to the distribution of the differences between the residuals, and can be min- 
imized by requiring that the slopes be small, as is done in the hit resolution measurement |7|. 
Any quantitative assessment is nevertheless difficult since it is not easy to disentangle the con- 
tributions of the rotations around the different axes. 

Only events with a single reconstructed track were used in this analysis. Module pairs with 
more than 35 (100) hits in the overlap region were analyzed in the pixel (strip) detectors. Only 
barrel layers satisfied these requirements. The lower hit requirement in the pixel modules 
reflects their smaller area and leads to correspondingly less precise mean values. 
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Figure 8: Relative shift between module pairs in the local u (top left) and local v (top right) 
direction in the BPIX, in the local w-coordinate in the TIB (bottom left), and TOB (bottom right). 
Only modules in the slice 80° < <p < 100° are shown in the TIB and TOB. 



The relative shifts are shown in Fig. [HJ In the strip detectors, the relative shift can be measured 
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Table 2: RMS of the mean of the distributions of the relative shift between overlapping module 
pairs, scaled by 1 / y/2 to account for the two independent measurements. Results are shown 
for the three data and two MC geometries described in the text. 
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in the local w-coordinate. In the pixel detector, this shift can be measured both in the local u and 
local v coordinates. The RMS of the distribution of the mean values are given in Table |2j where 
the values are scaled by 1/ a/2 to account for two independent measurements. The correspond- 
ing mean values of the distributions after alignment are 1.5 f/m or less. The expectation from 
MC, for both aligned as well as ideal geometry, is in good agreement with the distribution of 
the median of residuals quoted in Table [TJ where results in the pixel detector are more affected 
by the small size of the overlap sample. The differences in data may be a sign of yet unquan- 
tified systematic effects, like the aplanar distortions of modules. While the distributions of hit 
residuals taken over the full active region of the modules tend to be well centered, this is not 
always the case for distributions in the overlap regions, which are located at the edges of the 
active areas in the modules. However, this effect is smaller then 10 }im in the most precise u 
direction and is 13 ]dm in the V direction in the pixel barrel modules. There is also visible im- 
provement of the results with survey geometry compared to the geometry before alignment in 
the TIB, while in the TOB no module-level survey was done and in the BPIX no survey data 
between overlapping ladders were used either. 

5.4 Track parameter resolution 

The track parameter resolutions can be validated with independent reconstruction of the upper 
and the lower portions of cosmic ray tracks and comparison of track parameters at the point 
of closest approach to the nominal beamline. Both the upper and lower track segments were 
required to have at least three pixel hits. This mimics the topology of collision tracks. The track 
segments were reconstructed independently. Figure [9]shows the difference between upper and 
lower portions of tracks for all five track parameters. There is significant improvement due 
to tracker alignment, with good agreement between data and MC simulations. The results of 
the combined method are approaching those of a MC simulation with ideal detector geometry. 
The Pt measurements are most sensitive to the strip part of the tracker, while the other four 
parameters are dominated by the alignment of the pixel detector. The normalized distributions 
in Fig. [9] also show that the error estimates on the track parameters are in good agreement with 



predictions from MC simulations. Figure 10 shows the dependence of the resolution on pj for 
the track parameters d xy and 1/pr- 



5.4 Track parameter resolution 
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Figure 9: Differences between upper and lower track segment parameters measured at the 
point of closest approach to the beamline and scaled by l/v2- Distributions are shown for the 
distance of closest approach in the transverse direction d xy (top left), the same in the longitu- 
dinal direction d z (top right), the track azimuthal angle <p (middle left), the track polar angle 9 
(middle right), and 1/pj (bottom left). The plot on the bottom right shows the 1/pr difference 

normalized to its error, that is (1/ pj,i — l/pr,2)/< 



r l/p T ,l + a l/p T ,2 



Results are shown for four 



geometries: data before alignment (black dotted lines), data with combined method alignment 
(red solid), combined method MC (green dashed), and ideal MC (blue dash-dotted). 
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Figure 10: Dependence on pj of the differences between the track parameters measured at the 
point of closest approach to the nominal beamline, in the two halves of a cosmic ray track and 
scaled by l/y/2. The RMS of the distribution truncated at 95% is quoted at each momentum 
interval, shown for the distance of closest approach in the transverse direction d xy (left) and for 
1/fT (right). Results are shown for the geometries derived from ideal MC (blue squares), and 
from the alignment result with cosmic ray data (red circles) and combined method MC (green 
triangles). 



5.5 Systematic misalignment studies 

A global translation and rotation of the whole tracker is an example of a trivial transforma- 
tion which leaves the x 2 value of Eq. ([lj unchanged. This transformation has no effect on the 
internal alignment and is easily resolved by a suitable convention when defining the global 
reference frame. The convention that the center-of-gravity of all modules coincides with the 
design position was used in this work. A similar convention was used for the rotation of the 
tracker. 

There are, however, non-trivial transformations, so-called "weak" modes, that also preserve 
the x 2 value of Eq. ([lj and are of larger concern. The presence of weak modes in the geome- 
try resulting from track-based alignment was investigated following the approach described in 
Ref. |fl5|. Nine systematic distortions, in Ar, A(p, and Az as a function of r, (p, and z were applied 
to the aligned geometry. Studies were performed with both global and local alignment meth- 
ods and the results presented here were obtained with the global method. The systematically 
misaligned geometries were used as a starting point and the alignment procedure described in 
Section[4]was repeated. The nine geometries obtained after the alignments were then compared 
to the original aligned geometry to see if the distortions can be recovered by the alignment pro- 
cedure. 

The results were analyzed separately for the pixel (BPIX and FPIX), barrel strip (TIB and TOB), 
and forward strip (TID and TEC) sub-detectors. The remaining displacements of the modules 
in the TIB and TOB, which are expected to have the best illumination from cosmic ray tracks, 
are shown in Fig. [II] for four systematic deformations: the layer rotation (Acp = C\ + c^r), the 
twist deformation (A(p = c\z), the skew (Az = C\ cos <p), and the z-expansion (Az = C\z). The 
layer rotation can be well recovered in alignment with cosmic ray tracks, while the scatter in 
the module position difference is an evidence for weak modes in other projections. The twist 
and the skew deformations are reflected to a lesser degree in track x 2 , and therefore are only 
partially recovered. Finally, as the distribution of the track x 2 changes only marginally with 
respect to the z-expansion, this deformation cannot be recovered using cosmic ray tracks only. 



5.5 Systematic misalignment studies 
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Figure 11: Comparison of the position of the TIB and TOB modules with respect to the geom- 
etry obtained with the global method after applying systematic distortions (black solid lines) 
and after alignment (red dots) is shown on the left for layer rotation (top row), twist (second 
row), skew (third row), and z-expansion (bottom row) weak modes. For the twist misalign- 
ment, the shorter and longer solid lines correspond to the modules in the inner layer of the TIB 
and in the outer layer of the TOB, respectively. The plots on the right show the distributions 
of the corresponding track x 2 1 ndf after the alignment with the global method (blue solid line), 
after introducing the systematic misalignment (black dashed line), and after re-aligning (red 
solid line, below the blue solid line). 
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Radial deformations are generally recovered using constraints from overlapping modules. Also, 
mid-plane differences between the upper and lower legs of cosmic ray tracks are sensitive to 
deformations that differently affect the upper and lower parts of the tracker, including the 
layer rotation and telescope (Az = C\r) modes. However, more subtle deformations, including 
those discussed above, may be difficult to recover with cosmic ray tracks alone. Additional 
information with tracks from LHC beam interactions, which should be uniform in (p, and with 
additional constraints from the verticies and the masses of resonances, may provide better sen- 
sitivity to those systematic deformations. In principle, the Laser Aligment System (Section |5. 7} 
and hierarchical survey measurements provide complementary information, but evaluating 
systematic biases in those measurements is a challenging task. 

5.6 Study of detector geometry 

Two sets of comparisons between geometries obtained from track-based alignment with data 
have been performed: a comparison between the geometries from the local and the global 
methods, and a comparison between the geometry obtained with the combined method with 
respect to the design one. A comparison between two geometries is done after correcting for 
an overall residual shift and rotation of the whole detector, or a sub-detector, with respect to its 
center-of-gravity. 



Figure 12 shows the distribution of the differences between the rep positions of the 768 BPIX 
modules obtained with the local and the global methods. The RMS of the distribution is about 
12 ^m. The distribution of the differences between the z positions is similar once the two main 
weak modes, Az vs. z and Az vs. </>, are properly identified and removed by the means of the fit 



with the functions shown in the left of Fig. 11 Similar module-by-module comparisons were 
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performed on larger structures, like the TIB or TOB, where weak modes cannot be easily re- 
moved, and on coordinates for which the track-based alignment with cosmic rays has a limited 
sensitivity. This led to distributions characterized by a larger spread. This is also confirmed by 
alignment on simulated data and limits the accuracy of the absolute position determination. 

Assuming that the geometry of the tracker is best described by the one obtained with the com- 
bined method, a comparison with the design geometry indicates that the two BPIX half-barrels 
are shifted along the vertical axis by about 0.4 mm and the two half-barrels of the TIB have an 
extra separation along the z axis of about 5 mm. Both displacements are mechanically possible. 



5.7 Validation with the Laser Alignment System 
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The large displacement of the TIB half-barrels is supported by the optical survey measurements 



described in Ref. [9j. Figure 13 shows the displacement in the y- and in the z-coordinates of the 
BPIX and TIB modules with respect to the design position. In the case of the TIB, a large scatter 
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Figure 13: Comparison of the position of the modules in the combined method geometry with 
respect to the design one: {y design — ycomb) for the BPIX modules (left) and {zd es i sn — z CO mb) for 
the TIB modules (right) as a function of (p. In the plot on the left crosses (circles) represent 
modules of the positive (negative) x BPIX half-barrel. In the plot on the right crosses (circles) 
represent modules of the positive (negative) z TIB half-barrel. 



of the z-coordinate of the modules with respect to their design position is observed. The mod- 
ulation of the displacement as a function of the cp position of the module can be explained by 
the presence of the weak mode skew in the combined geometry as discussed in the previous 
section. 

5.7 Validation with the Laser Alignment System 

An independent test of the silicon tracker alignment was provided by the Laser Alignment 
System (LAS), which uses a system of 40 infrared laser beams (A = 1075 nm) to survey the po- 
sition of the large-scale structure elements of the tracker. The LAS measurements are available 
for 434 silicon strip modules, which are distributed over eight azimuthal sectors. The light is 
detected directly on the active area of the silicon sensors and therefore provides excellent beam 
position resolution with respect to the module sensitive area. For each TEC disk, there are 16 
modules distributed uniformly in cp at two radial positions that are intersected by laser beams. 
Furthermore, the LAS is capable of measuring the relative orientation of both endcaps and the 
TIB and TOB half-barrels. The TID and the pixel detector are not included in the LAS. 

Operation of the LAS during data-taking was mainly devoted to commissioning. In total, 31 out 
of the 40 laser beams were fully operational during this period. The available data were used 
to calculate the relative positions of the TEC disks. Using LAS measurements in the alignment 
procedure by means of additional terms to Eq. ^ is not discussed in this paper. 

For the calculation of the alignment parameters from LAS data, the TEC disks were treated as 
solid objects, each with one rotational and two translational degrees of freedom. Additional 
degrees of freedom resulted from the fact that the beam directions may deviate from their 
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nominal values at the level of 1 mrad. The calculation was based on an analytical approach for 
solving the ^-minimization problem for the laser hit residuals. 

For comparison, the position of the TEC disks was also calculated from the results of the track- 
based alignment using the combined method, where only the position of the modules illumi- 



nated by the LAS beams were considered, for consistency. Figure 14 compares the results of the 
measurements for positive z side TEC disks from the cosmic ray data analysis, the LAS data 
analysis, the LAS sector test [9|, and optical survey. All data are in good agreement. However, 
it is difficult to estimate the systematic uncertainties in the results for the track-based alignment 



due to the effects discussed in Section 5.5 Small movements of particular disks are expected 
between the sector test and optical survey measurements and the CRAFT results from the LAS 
and the track-based alignment. The optical surveys were done with the endcaps in a vertical 
orientation whereas they have a horizontal orientation within the tracker. 
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Figure 14: The three alignment parameters for positive z TEC disks measured using cosmic 
ray data (solid blue circles), by the LAS in this analysis (open black circles), during the LAS 
sector test [9] (red squares), and by optical survey (green triangles): rotation around global z 
(left) and translation in global x (middle) and in global y (right). 



As a comparison to track-based results, the laser hit residuals, defined as the difference between 
the nominal and the measured laser hit position along the ^-coordinate, have been determined 
from the tracker geometry before alignment as well as assuming the geometry from the com- 
bined method. Corrections have been made for deviations of the beams from their nominal 
orientations. The results are shown in Fig. 15 Assuming the aligned geometry as measured 
with tracks, the overall RMS of the residual distribution decreases significantly with respect 
to the before alignment case (cf. Table |3|. For the modules in the TEC, the RMS decreases by 
more than a factor of two, while the effect in the TOB is slightly smaller. The largest relative 
improvement is observed for the TIB modules. However, the resulting RMS of 200 ^rad for 
TIB modules remains larger than in the other sub-detectors. The observed residual widths may 
reflect an insufficient calibration of the lasers and a non-optimized laser beam correction pro- 
cedure. Furthermore, weak modes for which the track-based alignment with cosmic rays has 
limited sensitivity may contribute to the residual distributions observed in the LAS studies. 



6 Summary and Discussion 

In summary, results of the first full alignment of the CMS tracker have been presented. Two 
algorithms have been used to determine the positions of all 16 588 silicon modules. The two 
alignment methods have been combined sequentially to take into account most effectively both 
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Table 3: RMS of laser hit residuals (^nominal — ^measured) for the TEC, TIB, and TOB before and 
after alignment with the combined alignment method. 
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Figure 15: Laser hit residuals ((/'nominal — ^measured) for the TEC (left), TIB (middle), and TOB 
(right), using the geometry before alignment (dashed black line) and the geometry obtained 
from the combined alignment method (solid red line). 



global and local correlations of module positions. The results are based on the analysis of about 
three million cosmic ray tracks recorded with a 3.8 T magnetic field. 

The precision of the detector positions with respect to particle trajectories has been derived 
from the distribution of the median of the cosmic muon track residuals to be on average 3- 
4 ^m RMS in the barrel and 3-14 RMS in the endcaps in the most sensitive coordinate. 
These results are supported by the difference of the residuals in regions where modules overlap 
within a layer, thus reducing effects of multiple scattering in the calculation of the residuals. 
Nevertheless, additional small systematic effects may still be present. The achieved tracker 
resolution in all five track parameters has been checked with a study of the two halves of 
a cosmic ray track and compared with predictions from a detailed detector simulation. The 
measured resolutions are close to those that would be observed in a detector with perfectly 
placed modules. 

Overall, a significant improvement in track fit performance compared either to using the be- 
fore alignment or survey geometry has been observed. The resulting aligned geometry allows 
a study of the assembly precision of the individual sub-detectors. A comparison of module 
positions obtained with complementary methods supports findings based on track residuals. 
Nonetheless, certain deformations in the geometry that do not change the % 2 of the tracks, 
cannot be ruled out due to the largely vertical nature of the cosmic track data. 

Clear improvement in the LAS residuals is observed when using the result of the track based 
alignment. The stand-alone LAS data analysis agrees well with that from previous tests and 
with survey alignment. The operation of the LAS has shown that the laser beams operate 
properly. 
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Experience gained in the alignment analysis of the silicon modules with cosmic ray particles 
is valuable in preparation for the CMS tracker alignment with the data from LHC collisions, 
which is critical to achieving the physics goals of the CMS detector. Integration of measure- 
ments from cosmic ray and collision tracks, LAS, and survey will be critical for the optimal 
tracker alignment with the first data expected from LHC beam collisions. Because the track pa- 
rameter resolutions are close to the design performance with cosmic ray data already, statistical 
uncertainties in track parameters are not expected to be considerably affected by alignment at 
CMS startup with data from beam collisions. 
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